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1  INTRODUCTION 

1.1.  Research  Background 

Heat  generated  by  electronic  devices,  laser  arrays,  and  power  sources  in  Missile  Defense  Agency  (MDA) 
facilities  needs  to  be  dissipated  effectively  to  ensure  that  temperature-sensitive  components  function 
properly  and  achieve  their  maximum  performance.  However,  the  thermal  management  needs  at  MDA  are 
often  more  challenging  than  those  in  industry  because  of  weight  and  volume  limitations  in  the  platforms. 
To  remove  heat  dissipation  at  the  rate  above  10^  W/m^,  the  heat  transfer  requirement  is  beyond  the 
capabilities  of  the  conventional  cooling  methods.  MicroChannel  (1  pm  <  hydraulic  diameter  <  1  mm) 
heat  sink,  introduced  by  Tuckerman  and  Pease^^^  in  the  early  1980s,  is  an  innovative  cooling  technology, 
which  could  achieve  very  high  rate  of  energy  flow  by  forcing  the  cooling  fluid  through  microchannels 
fabricated  directly  in  the  silicon  wafer  containing  heat-dissipating  large  scale  integrated  circuit  (LSIC) 
components.  This  new  technology  demonstrates  two  advantages  over  the  conventional  heat  transfer 
methods.  The  first  one  is  that  microchannel  can  be  manufactured  near  where  the  heat  source  is  located 
without  impacting  the  mechanical  structure  and  therefore  is  capable  of  catching  heat  where  it  is 
generated.  The  shorter  heat  conduction  path  between  the  heat  source  and  the  sink  in  this  microchannel 
heat  sink  configuration  results  in  a  much  lower  thermal  resistance.  The  other  advantage  of  the 
microchannel  heat  sink  is  its  superior  cooling  capacity  due  to  the  great  surface  area  per  unit  volume  of 
microchannels. 

Even  though  microchannel  heat  sinks  are  capable  of  dissipating  much  higher  heat  fluxes  than 
conventional  cooling  techniques,  two  existing  problems  are  limiting  their  applications.  The  first  one  is  the 
large  pressure  drop  along  the  flow  direction  caused  by  the  microscale  hydraulic  diameter,  which  requires 
more  pumping  power.  The  second  problem  is  that  the  low  flow  rate  in  microchannel  produces  large 
temperature  rise  along  the  flow  direction  in  both  the  cooling  fluid  and  the  solid  walls,  which  may  damage 
temperature-sensitive  electronic  components  of  LSIC.  The  problem  becomes  worse  if  temperature 
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distribution  on  a  plate  or  a  printed  circuit  board  is  highly  biased. 

To  fully  utilize  the  advantages  of  microchannel  cooling  technology,  such  as  high  energy  mitigation 
capability  and  low  heat  resistance  and  effectively  solve  the  existing  two  problems  encountered  by  plane 
microchannels,  i.e.,  high  pressure  drop  and  high  channel-wise  temperature  rise,  this  research  project  is 
going  to  apply  vortex  generator  (VG)  technology^^'"^^  in  microchannels,  i.e.,  periodically  mounting  arrays 
of  square  cylinders/bars  in  microchannel  to  generate  self-sustained  oscillating  flow.  Since  the  ordered 
laminar  oscillating  flow  requires  much  less  pumping  power  than  turbulent  flow  to  achieve  the  same  heat 
transfer  rate  or  to  get  a  much  higher  heat  transfer  rate  at  the  same  pumping  power,  due  to  the  advantage 
of  ordered  laminar  self-sustained  oscillation  over  random  turbulent  fluctuations  which  yields  less  viscous 
dissipation.  Meanwhile,  the  strong  convective  mixing  among  the  communicating  channels  helps 
reducing  the  channel-wise  temperature  rise.  A  numerical  simulation  was  reported  that  this  VG  equipped 
microchannel  yields  much  higher  cooling  capacity  than  plane  microchannels  without  causing  a  large 
pressure  drop^^^.  This  new  technology  is  being  anticipated  to  be  used  in  MDA  applications  to  design 
more  efficient  and  more  reliable  microchannel  heat  sinks  or  microchannel  heat  exchangers  with  high  heat 
dissipation  capabilities  (>10^  W/m^)  and  low  heat  resistance  in  order  to  meet  the  increasing  heat 
dissipation  demand  for  MDA  electronic  devices  at  ground,  air,  and  space  platforms  at  both  ambient/room 
and  cryogenic  temperatures. 

To  the  best  of  the  investigators’  knowledge,  currently  there  are  no  available  experimental  data  in  the  open 
literature  that  could  lead  to  a  practical  design  for  self-sustained  oscillating  flow  in  microchannels. 
Furthermore,  traditional  numerical  approaches,  such  as  those  using  Navier-Stokes  equations,  are  not  able 
to  correctly  predict  the  hydraulic  and  thermal  behaviors  of  microchannel  flows.  The  experimental  results 
have  shown  that  fluid  flow  and  heat  transfer  characteristics  of  microchannels  deviate  from  that  of 
conventionally  sized  channels  (macroscale,  hydraulic  diameter  >  1  mm)^^'^l  For  liquid  flow  in 
microchannels,  the  flow  characteristics  are  more  complicated  than  that  of  gas.  In  addition  to  the 
rarefaction  effect,  effects  of  electric  double  layer  (EDL)^^^^  and  dissolved  noncondensable  vapor^^^^  could 
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account  for  a  certain  percentage  (about  the  order  of  10%)  of  the  heat  transfer  rate  variations  in 
microchannel  flows,  which  were  usually  negligible  in  macroscale  channel  flows.  To  date,  there  have 
been  no  acceptable  models  for  simulating  liquid  flow  in  microchannels  that  take  into  accounts  all  of  the 
effects  due  to  the  microscale  size. 

The  available  investigations  on  the  physical  mechanism  of  heat  transfer  enhancement  by  oscillating  flow 
was  primarily  concentrated  on  VG  equipped  conventionally  sized  channels  through  both  experimental 
and  numerical  methods  It  was  found  that  the  heat  transfer  rate  and  pressure  drop  are  both 

increasing.  The  physical  mechanism  of  heat  transfer  enhancement  for  VG  equipped  conventionally  sized 
channels  was  concluded  as  follows 

•  Surface  interruption  prevents  the  continuous  growth  of  the  thermal  boundary  layer  by  periodically 
interrupting  it.  This  mechanism  occurs  even  at  low  Reynolds  numbers  when  the  flow  is  steady  and 
laminar; 

•  Above  a  critical  Reynolds  number,  the  interrupted  surfaces  offer  an  additional  mechanism  of  heat 
transfer  enhancement  by  inducing  oscillations  in  the  flow  in  the  form  of  shed  vortices.  These 
vortices  enhance  local  heat  transfer  by  continuously  bringing  fresh  fluid  towards  the  heat  transfer 
surfaces; 

•  In  addition  to  heat  transfer  enhancement,  surface  interruption  increases  the  friction  loss  and  thus 
requires  larger  pumping  power.  This  is  partly  due  to  the  larger  skin  friction  associated  with  the 
hydrodynamic  boundary  layer  restarting.  In  the  unsteady  flow  regime,  the  time-dependent  flow 
behavior  associated  with  vortex  shedding  increases  friction  loss  through  the  stokes  layer  dissipation 
and  form  drag  through  Reynolds  stresses; 

•  The  boundary  layer  restart  and  self-sustained  oscillatory  mechanisms  simultaneously  influence  both 
the  overall  heat  transfer  and  the  pumping  power  requirement. 
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While  for  the  investigation  on  VG  equipped  microchannel  flows,  only  one  numerical  study  is  available^^^. 
The  numerical  results  indicated  that  this  VG  equipped  microchannel  yields  much  higher  heat  transfer  rate 
without  causing  a  large  pressure  drop,  which  is  different  from  that  of  conventionally  sized  channels.  This 
numerical  observation  greatly  encourages  this  project  to  investigate  the  physical  mechanism  of  heat 
transfer  performance  and  flow  resistance  on  VG  equipped  microchannel  flows. 

In  this  project,  a  series  of  experimental  studies,  including  micro-particle  image  velocimetry  (PIV), 
hydraulic  and  thermal  measurements  will  be  carried  out  to  reveal  the  physical  mechanism  of  self- 
sustained  oscillating  flow  on  heat  transfer  and  flow  resistance  in  microchannels.  In  addition, 
comprehensive  parameter  studies  will  be  conducted  to  obtain  correlations  among  the  critical  control 
parameters,  such  as  the  Reynolds  number,  Nusselt  number  and  channel  geometric  parameters,  which  will 
become  unique  criteria  for  the  future  design  of  microchannel  heat  sinks. 

1.2.  Project  Objectives 

Self-sustained  oscillating  flow  in  microchannels  is  a  promising  new  concept  in  thermal  management. 
Once  it  is  well  understood,  it  can  be  used  for  the  cooling  of  electronic  devices  or  other  equipment  that 
requires  large  heat  transfer  flux  concerned  with  high  temperature  sensitive  components  in  MDA  related 
systems.  The  primary  objectives  of  this  project  are: 

•  To  understand  the  hydraulic  and  thermal  characteristics  of  self- sustained  oscillating  flow  through 
microchannels.  Flow  visualization,  flow  spectrum  study,  and  pressure  drop  and  temperature 
measurements  will  be  conducted  for  single-row  and  multi-row  channels  to  investigate  the  physical 
mechanisms  of  flow  and  heat  transfer.  For  the  single-row  channel,  the  focus  will  be  the 
characteristics  of  heat  transfer  and  pressure  drop.  For  the  multi-row  channel,  the  fluid  mixing 
processes  among  communicating  channels  will  be  the  research  focus  in  addition  to  the  heat  transfer 
and  flow  resistance  characteristics. 
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•  To  obtain  correlations  and  recommend  parameters  for  the  optimum  design  of  a  highly  efficient  heat 
sink  or  heat  exchanger.  Various  parameters,  such  as  streamwise  and  transverse-wise  pitches,  block 
ratio,  Reynolds  number  (and/or  Prandtl  number)  will  be  varied  to  investigate  their  influence  on  the 
target  parameters,  such  as  the  Nusselt  number  and  the  apparent  friction  coefficient.  The  generated 
database  will  become  a  unique  reference  for  the  development  or  design  of  a  highly  efficient  heat  sink 
or  heat  exchanger. 

1.3.  Significance  of  the  Project 

This  project  will  contribute  to  the  MDA’s  Innovative  Science  and  Technology  Program  goal  of 
establishing  advanced  thermal  management  approaches  to  enable  electronic  devices  to  achieve  their 
maximum  performance  with  reduced  pumping  power  and  temperature  sensitivity,  and  meanwhile  greatly 
enhance  Florida  International  University  (FIU)’s  research  and  development  capabilities  related  to  the 
needs  of  MDA,  which  will  lead  to  FIU’s  future  participation  in  the  primary  MDA  research  programs. 

With  an  enrollment  of  over  33,000  students,  including  about  70%  minority  students,  FIU  is  the  largest 
producer  of  Hispanic  engineering  graduates  in  the  nation.  It  provides  a  good  chance  to  expose  minority 
students  to  new  technology  development  in  the  field  of  thermal  management. 
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2  Experimental  Study 

2.1  Experimental  setup 

Figure  1  shows  the  picture  of  the  hydraulic  loop  used  for  the  experimental  investigations  of  this  project. 
Schematic  diagram  of  the  loop  is  shown  in  Figure  2  for  the  convenience  of  recognizing  each  component. 
In  the  experimental  study,  DI  water  was  pumped  through  a  5 -pm  filter  to  prevent  any  particles  from 
entering  the  microchannels  and  creating  blockages.  A  micro  flow  meter  was  installed  downstream  of  the 
filter  to  measure  the  flow  rate,  local  temperature  and  hydraulic  loop  pressure  of  the  flowing  water. 
Before  entering  the  test  section,  the  DI  water  was  directed  into  a  temperature  controlled  water  bath  to 
keep  the  inlet  water  temperature  at  a  constant  value  in  the  range  of  20-100^C  for  each  run  of  test.  The 
water  temperature  was  increased  after  flowing  through  the  heated  test  section  and  than  cooled  down  to  a 
certain  value  around  the  inlet  temperature  through  a  chill  water  heat  exchanger  located  downstream  of  the 
test  section  exit.  After  leaving  the  heat  exchanger,  the  DI  water  returned  to  the  pump,  completing  the 
hydraulic  loop. 
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Fig.  1  Photograph  of  the  hydraulic  loop. 
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Fig.  2  Schematic  diagram  of  the  hydraulic  loop. 
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Test  section 

Each  test  section  was  composed  of  three  plates/chips:  the  capping/bottom  plate  was  made  of  quartz  for 
visualization  study,  and  the  middle  and  top  plates  were  made  of  silicon  wafer.  VG  cube  arrays  were 
fabricated  in  the  middle  one  and  the  top  plate  was  used  for  mechanical  support  and  thermal  couple 
installations. 


Fig.  3  Cross-section  view  of  the  test  section. 

The  three  chips  were  bonded  together  with  epoxy  resin  to  form  a  flow  channel  (Figure  3)  and  connected 
with  two  plenums  for  the  convenience  of  temperature  and  pressure  drop  measurements  (Figure  4).  Two 
thermocouples  were  inserted  into  the  plenums  to  measure  the  inlet  and  outlet  temperatures,  respectively. 
The  test  section  pressure  drop  was  measured  by  the  pressure  transducer.  Five  pairs  of  thermocouples 
were  embedded  in  five  measurement  slots  fabricated  in  the  supporting  plate  to  measure  the  surface 
temperature  of  the  flow  channel.  The  empty  of  the  slots  was  filled  with  aluminum  nitride  powder  to 
maintain  constant  thermal  boundary  conditions.  The  thermocouples  were  covered  by  an  electrical 
insulating  layer  to  prevent  any  electrical  interference  with  the  heater.  The  heat  input  was  delivered 
through  an  electrical  heater  (0-1200  W)  attached  to  the  supporting  plate,  and  adjusted  through  a  variable 
transformer.  The  test  section  and  heater  were  well  insulated  to  minimize  the  heat  loss  to  the 
surroundings,  and  sandwiched  with  stainless  steel  plates  for  the  safe  assembling  with  the  plenums. 
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Connection  with  plenums 

Close  up  of  the  test  section  and  plenum  connection  is  shown  in  Figure  4.  The  assembled  test  section  was 
first  glued  horizontally  to  two  connecting  plates  by  inserting  the  test  section  into  the  slots  fabricated  in 
each  plate.  The  ending  of  the  slots  was  made  exactly  identical  to  the  height  and  width  of  the  test  section 
and  the  opening  was  cut  slightly  larger,  for  the  purpose  of  preventing  any  overflow  of  the  glue  through 
the  slots  into  the  micro-sized  flow  channels.  The  two  connecting  plates  were  fastened  tightly  to  the 
plenums  with  plastic  gaskets  to  make  the  connection  watertight.  The  sandwiched  test  section  was  rigidly 
supported  with  stainless  steel  plates  and  bolts  to  avoid  any  stresses  or  loads  applied  on  the  microchannel 
chips. 


Fig.  4  Test  section  and  plenum  connection. 


Verification  of  Instrumentation  Calibration 

After  the  loop  was  checked  for  leakage,  calibration  of  individual  devices  and  components  was 
tested  against  the  calibration  results  provided  by  the  vendors. 
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Differential  pressure  transducer  was  tested  by  measuring  the  pressure  drop  in  a  straight  tube  (33”  long) 
oriented  horizontally  and  comparing  measurements  with  the  calculated  results  from  the  Darcy-  Weisbach 
equation: 


AP  =  f  — 


2 


(1) 


Where  /  is  friction  factor  and  Dh  is  the  hydraulic  diameter. 

The  comparison  results  are  shown  in  Figure  5.  It  can  be  seen  that  the  pressure  transducer  readings  are 
larger  than  the  theoretical  values  at  large  flow  rates,  which  could  be  due  to  the  tube  flow  being  not  fully 
developed  in  the  whole  measurement  length,  while  the  friction  factor  is  smaller  for  fully  developed  flow 
which  is  assumed  in  the  theoretical  calculations.  The  differential  pressure  transducer  is  calibrated  by  the 
vendor  to  an  accuracy  of  ±0.25  %  of  the  full  measurement  range  (Figure  6),  which  is  not  repeatable  or 
attainable  with  the  present  experimental  devices  and  is  assumed  to  be  acceptable. 


Figure  5.  Pressure  drop  comparison 
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Figure  6.  Vendor  provided  calibration  results  for  the  differential  pressure  transducer. 

The  flow  meter  was  tested  by  measuring  the  volume  of  water  with  a  measuring  jar  and  the  time  with  a 
stop  watch.  The  measured  results  are  compared  with  the  flow  meter  readings  in  Figure  7.  The  deviation 
of  the  readings  is  up  to  22%  of  the  measured  value  at  low  flow  rate,  and  within  ±5%  at  larger  one.  This 
result  is  close  to  the  accuracy  of  ±  2%  provided  by  the  vendor. 


Figure  7.  Flow  meter  tests 
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All  thermocouples  have  been  tested  at  two  different  temperatures,  0  ^C,  provided  by  ice-water  mixture 
and  100  ^C,  provided  by  boiling  water,  at  atmospheric  pressure  (Figure  8).  The  accuracy  is  within  ±0.2 
^C,  which  is  the  calibration  accuracy  provided  by  the  vendor. 


Figure  7.  Thermocouple  calibrations 


2.2  Data  collection  and  reduction 


A  Lab  VIEW  data  acquisition  system  was  used  for  data  collection.  After  the  loop  has  been  running  for  a 
certain  period  of  time,  steady  state  was  assumed  when  the  temperature  variation  was  less  than  ±1% 
within  5  minutes.  All  temperature,  flow  rate  and  pressure  drop  data  were  then  recorded. 

The  hydraulic  diameter,  dh  of  the  test  section  is  defined  as: 


Where  Ac  is  the  cross-section  area  of  the  microchannel  and  P  is  the  perimeter. 
The  Reynolds  number.  Re  is  given  by: 


Re  = 


pud^ 


A 


(2) 


(3) 
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Where  u  is  the  mean  velocity  of  water;  fi  is  the  dynamic  viscosity  and  p  is  the  density.  The  reference 
temperature  for  properties  is  taken  as  the  average  of  the  wall  temperatures  and  the  fluid  inlet  and  outlet 
temperatures. 

The  Nusselt  number,  Nu  is  given  as: 


Nu  = 


hd^ 

k 


Where  h  is  the  heat  transfer  coefficient  of  water  and  k  is  the  thermal  conductivity. 
The  friction  factor, /is  calculated  from  the  pressure  drop  across  the  micro  channel: 


(4) 


/  = 


lAPd, 

Lpu^ 


Where  L  is  the  microchannel  length  and  AP  is  the  pressure  drop. 


(5) 


2.3  Experimental  results  and  discussions 


Upon  successful  testing  and  calibration  of  the  hydraulic  flow  loop  heat  transfer  studies  with  the 
VG-equipped  microchannels  were  conducted. 

Figures  8  and  9  present  the  variations  of  Nu  and  /  with  respect  to  Re  for  the  VG-equipped 
microchannel  with  one  row  of  40  VG  cubes.  Other  dimensions  of  this  microchannel  are  as  follows:  height 
505  pm,  length  39895  pm,  and  width  2525  pm.  Experimental  results  are  compared  on  the  same  figures 
with  the  correlations  developed  for  smooth  microchannels 
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C 

r  _  f  ,turb 


(9) 


Where  Pr  is  the  Prandtl  number;  Wc  is  the  mieroehannel  eenter  to  eenter  spaeing;  Z  is  the  dimensionless 
parameter  C/is  the  empirieal  eoeffieient.  Subseripts  lam  stands  for  laminar  and  turb  for  turbulent. 

From  Figures  8  and  9,  it  ean  be  observed  that  with  one  row  of  VG  eubes  mounted  in  the 
mieroehannels,  Nu  is  larger  than  that  predieted  by  laminar  eorrelation  and  smaller  that  predieted  by 
turbulent  one  at  7?e<630,  where  the  flow  regime  is  not  determined  in  the  present  study.  The  /  is  larger 
than  the  laminar  eorrelation  predietions  at  Re  >  400,  and  also  larger  than  the  turbulent  eorrelation 
predietions  at  Re  >  670. 


Figure  8.  Nu  vs.  Re  for  the  1x40  cube  array  mieroehannel. 
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As  shown  in  Figure  fRe  is  much  larger  and  oscillating  with  Re.  The  temperature  ratio,  shown  in 
Figure  11,  shows  that  the  temperature  increases  more  in  the  water  flow  direction  than  that  on  the  wall, 
indicating  the  improvement  of  wall  temperature  distribution  by  VG  cubes. 

These  experimental  results  indicate  that  VG  cubes  could  improve  the  heat  transfer  and  wall 
temperature  distribution  in  microchannel  compared  with  smooth  one,  while  the  flow  resistance  is  larger. 
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The  effects  of  the  inlet  fluid  temperature  and  the  heat  input  have  been  studied  for  a  multi-row 
VG-equipped  microchannel.  This  channel  has  7  x  21  array  of  VG  cubes,  its  other  dimensions  are  as 
follows:  height  505  pm,  length  39895  pm,  and  width  11615  pm.  Three  different  inlet  fluid  temperatures 
have  been  tested  for  two  different  heat  inputs.  Figure  12  presents  the  variations  of  Nu  with  respect  to  Re. 
As  can  be  seen  from  the  graphs,  with  the  increase  of  inlet  fluid  temperature  the  Nusselt  number  is 
increased  for  the  given  test  conditions.  It  should  be  noted  that  with  the  increase  of  Reynolds  number,  the 
Nusselt  number  doesn’t  vary  much.  It  can  also  be  noted  that  the  Nusselt  number  is  almost  the  same  for 
the  different  heat  input  conditions  with  the  same  inlet  fluid  temperature.  This  could  further  prove  that  the 
Nusselt  number  doesn’t  change  much  with  change  in  heat  flux. 

For  a  flxed  Re  =  200,  the  effects  of  change  in  inlet  fluid  temperature  on  the  Nusselt  number  is 
given  in  Figure  13.  There  was  an  average  increase  of  40  %  in  the  Nusselt  number  when  the  temperature 
of  the  inlet  fluid  has  been  increased  by  5  ®C. 
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Figure  12.  Nu  vs.  Re  for  the  7x21  cube  array  microchannel. 
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Figure  14  shows  the  heat  transfer  and  frietion  faetor  results  with  variation  of  input  power  to  the 
mieroehannel  with  a  eonfiguration  of  7  rows  and  21  eolumns  of  VG  bars.  To  study  the  effeet  of  test 
parameters,  the  data  has  been  statistieally  analyzed  and  it  has  been  found  that  the  data  is  normally 
distributed  and  the  heat  flux  and  Reynolds  number  and  heat  flux  variation  has  a  signifieant  effeet  on  the 
heat  transfer  rate  while  the  ehange  in  inlet  temperature  of  the  fluid  does  not  have  a  signifieant  effeet.  The 
Nusselt  number  ean  be  seen  to  inerease  signifieantly  after  Reynolds  number  of  100.  In  the  range  of 
Reynolds  number  tested,  an  inerease  in  10  watts  gave  about  15  %  inerease  at  Re  «  100  and  11  %  at  Re  « 
500  in  Nusselt  number. 
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Figure  14.  Nu  and / vs.  Re  for  the  7x21  cube  array  microchannel. 


In  Figure  15,  the  near-linear  effeet  of  heat  flux  level  on  Nusselt  number  ean  be  seen  for  a  flxed 
inlet  temperature  of  20®C.  The  Reynolds  number  was  kept  eonstant  for  eaeh  of  the  four  different  eases.  It 
was  observed  that  at  flxed  heat  flux,  the  Nusselt  number  ehanged  by  about  5%  with  ehange  in  Reynolds 
number  of  100  within  the  range  of  100-500. 


Figure  15.  Nu  vs.  q  for  the  7x21  cube  array  microchannel. 
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We  define  the  Colburn  j  factor  as  a  measure  of  heat  transfer 


Nu 

RePr®" 


(10) 


The  Colburn  j  factors  obtained  from  the  experiment  are  found  to  be  larger  than  those  obtained  by  steady 
and  unsteady  simulation  models  of  Zhang  et  al^^^^,  confirming  that  the  overall  heat  transfer  was  improved. 
Figure  16  shows  the  Colburn  j  factors  for  four  heat  fluxes  used  in  the  experiments. 


Figure  16.  Colburn  j  factor  comparison  for  the  7x21  cube  array  microchannel. 


In  Figure  17,  the  friction  factor  is  plotted  against  the  Reynolds  number  for  the  7  x  21  cube  array 
microchannel.  The  friction  factor  for  all  the  cases  is  seen  to  follow  a  certain  trend.  The  product  /Re  and 
temperature  ratio  response  with  change  in  Reynolds  number  are  shown  in  Figures  18  and  19. 
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Figure  17.  /  vs.  Re  for  the  7x21  cube  array  microchannel. 
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Figure  18.  fRe  vs.  Re  for  the  7x21  cube  array  microchannel. 
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Figure  19.  Temperature  ratio  vs.  Re  the  7  x  21  cube  array  microchannel. 

Figure  20  shows  the  surfaee  temperature  variation  along  the  length  of  the  ehannel.  It  ean  be  seen 
that  at  high  Re  the  surfaee  temperature  is  deereased,  whieh  suggests  more  heat  transfer.  But  the  Nusselt 
number  variation  is  not  mueh.  This  may  be  attributed  to  the  faet  that  the  heat  transfer  is  not  only  taking 
plaee  by  eonveetion  but  also  by  eonduetion.  Thus,  eonduetion  as  well  as  eonveetion  heat  transfers  may  be 
inereasing  by  similar  proportions  when  the  Reynolds  number  is  inereased. 

Figure  21  shows  the  inlet,  outlet  and  wall  temperature  distributions  along  the  test  seetion  with  7  x 
21  eube  array.  The  inlet  temperature  was  kept  eonstant  at  20  and  the  heat  input  was  maintained  at  7 
Watts.  Results  at  Re  of  200  and  300  are  given.  It  ean  be  seen  that  the  wall  temperature  was  deereased 
along  flow  direetion  at  higher  Reynolds  number.  While  the  ehange  in  Reynolds  number  had  an 
insignifieant  effeet  at  the  entranee  of  the  test  seetion,  there  was  5  %  deerease  in  the  wall  temperature  at  a 
distanee  of  40  mm  from  the  entranee. 
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Figure  20.  Variation  of  surface  temperature  along  the  7x21  cube  array  microchannel. 
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Figure  21.  Variation  of  surface  temperature  along  the  7x21  cube  array  microchannel. 
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Figure  22  shows  the  temperatures  ehange  with  time  at  a  distanee  of  15  mm  (a)  and  40  mm  (b)  for 
Reynolds  number  at  200  and  500.  Attempts  were  made  to  verify  if  there  is  an  oseillation  of  the  wall 
temperatures.  The  temperature  was  reeorded  for  every  100  milliseeonds.  There  were  some  fluetuations 
observed  and  whieh  were  of  a  magnitude  of  ±  0.5%  from  the  average  value.  Most  of  them  were  found  to 
be  random  fluetuations  rather  than  periodie. 


Figure  22.  Variation  of  temperature  with  time. 


2.4  Conclusions  from  the  experimental  study 

Heat  transfer  study  with  single-  and  multi-row  VG-equipped  mieroehannels  resulted  in  following 
findings: 

For  a  single-raw  VG-equipped  mieroehannel  the  heat  transfer  and  wall  temperature  distribution  in 
mieroehannel  are  improved  eompared  with  a  ehannel  without  VG  bars,  while  the  flow  resistanee  is  larger. 
For  a  multi-raw  VG-equipped  mieroehannel: 

•  inerease  of  inlet  fluid  temperature  results  in  inereased  Nusselt  number 

•  the  Nusselt  number  is  almost  the  same  for  the  different  heat  input  eonditions  with  the  same  inlet 
fluid  temperature 
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•  the  Nusselt  number  inereases  signifieantly  when  the  Reynolds  number  inereases  above  100 

•  the  Nusselt  number  inereases  almost  linearly  with  inereasing  heat  flux  to  the  mieroehannel 

•  experimentally  obtained  Colburn  j  faetors  are  larger  than  those  obtained  by  steady  and  unsteady 
simulation  models  of  Zhang  et  al^^^^,  eonfirming  that  the  overall  heat  transfer  was  improved 

•  with  inereasing  Reynolds  number  the  mieroehannel  surfaee  temperature  is  deereasing 

•  the  wall  temperature  was  deereased  along  flow  direetion  at  higher  Reynolds  number,  while  the 
ehange  in  Reynolds  number  had  an  insignifieant  effeet  at  the  entranee  of  the  mieroehannel 

•  no  mieroehannel  wall  temperature  oseillations  were  observed  under  studied  range  of 
experimental  eonditions 
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3  Numerical  Study 

3.1  Mathematical  models  and  numerical  methods 

3.1.1  Computational  domain 

Three  different  arrays  of  VG  bars  in  the  same  ehannel  were  used  to  numerieally  study  the  water  flow  and 
heat  transfer  in  a  3-D  mieroehannel,  with  H  =  W=  800  pm,  Dh  =  800  pm,  and  the  ehannel  length  over  the 
hydraulie  diameter  ratio  of  L/Dk  =  5.  Figure  23  (a),  (b)  and  (e)  presents  the  ehannel  geometry  with  the  1,  2 
and  4  VG  bars,  respeetively.  The  upstream  and  downstream  points  are  defined  for  traeing  the  history 
value  of  flow  properties  with  time. 


Upstream  point  Downstream  point 
(a) 
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Up  stream  point  Downstream  point 


(C) 

Figure  23.  MicroChannel  geometry:  a)  1-bar  array;  b)  2-bar  array;  c)  4-bar  array. 


3.1.2  Governing  equations 

Three-dimensional  eontinuity,  momentum  and  energy  equations  in  Cartesian  tensor  notations  are 
deseribed  as  follows: 


27 


Contract  Number:  FA8650-04-C-2405 


Continuity  equation 
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Energy  equation 
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3.1.3  Boundary  conditions 

The  schematic  for  boundary  conditions  is  shown  in  Figure  24.  The  boundary  conditions  are  described  as 
follows: 

Inlet  and  outlet:  Periodic  conditions,  specifying  the  mass  flow  rate  {Q  =  constant)  and  the  temperature. 
Wall:  All  the  surfaces  are  treated  as  no-slip  boundaries.  And  the  walls  were  adiabatic  ^”=0,  except 
that  the  bottom  wall  is  heated  with  constant  heat  flux^”=  const . 

The  velocities,  pressure  and  temperature  in  the  streamwise  direction  of  the  computational  domain  is 
considered  for  periodic  boundary  conditions.  The  assumption  of  periodicity  implies  that  the  velocity 
components  repeat  themselves  in  space  as  follows: 

u^  (r )  =  u^  {r  +L)  =  u.  (r  +  2L)  =  •-  (15) 
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Side  walls  Top  (surface)  wall 


Figure  24.  Schematic  for  boundary  conditions 

The  pressure  is  not  periodic.  Instead,  the  pressure  drop  between  modules  is  periodic: 

A/)  =  /)(?)-/»(?  +  £)=  +  £)-/)(f  +  2Z)=  •••  (16) 

Similarly,  the  temperature  drop  between  modules  is  periodic: 

Ar  =  r(r)-r(r+Z)=r(r+Z)-r(r  +  2Z)=---  (17) 

Where  r  is  the  position  veetor  and  L  is  the  periodie  length  veetor  of  the  domain  eonsidered. 

3.1.4  Model  description 

The  time-dependent  governing  Navier-Stokes  equations  were  solved  with  finite  volume  method.  The 
eontinuity,  momentum  and  energy  equations  were  solved  by  seeond  order  upwind  implieit  seheme  using 
an  unsteady  laminar  solver,  and  the  pressure-veloeity  eoupling  was  dealt  with  by  SIMPLEC  algorithm. 
For  the  eonjugate  heat  transfer  problem,  the  same  energy  equation  is  used  for  the  fluid  domain  as  well  as 
the  solid  domain  with  high  eonduetivity  and  zero  veloeity  within  the  VG  bars.  To  avoid  the  divergenee, 
the  under-relaxation  teehnique  was  applied  in  all  dependent  variables.  In  the  investigation,  the  under¬ 
relation  faetor  for  the  pressure  was  0.7,  that  for  the  veloeity  eomponents  was  0.9,  and  that  for  energy 
eonponent  was  0.9.  In  the  segregated  solution  method,  eaeh  diserete  governing  equation  was  linearized 
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implicitly  with  respect  to  that  equation’s  dependent  variable.  The  convergence  criterion  for  residuals  was 
set  to  10'^  for  all  the  continuity,  momentum,  and  energy  equations. 

3.1.5  Mesh  independence  study 

In  order  to  achieve  a  grid  independent  solution  for  the  microchannel  simulations,  computations  have  been 
performed  at  three  different  meshes  16  x  16  x  80,  20  x  20  x  100,  and  32  x  32  x  160  at  Reynolds  number 
of  300  by  comparing  the  values  of  averaged  Nusselt  number  on  the  heated  wall  and  the/^7?e  number.  The 
Nusselt  number  shows  variations  of  9.9  %  and  4.1  %  while  the/^7?e  number  reveals  variations  of  4.0  % 
and  1.0  %  for  16  x  16  x  50  and  20  x  20  x  100  grids  respectively  with  respect  to  the  grid  size  of  32  x  32  x 
160.  The  Nusselt  number  variation  with  mesh  number  can  be  found  in  Figure  25  (a)  and  (b) 

respectively.  Consequently  the  mesh  size  of  20  x  20  x  100  was  selected  as  grid  independent  result  and 
was  used  for  all  the  simulations.  Figure  26  shows  the  mesh  used  in  the  simulations  with  the  domain  with 
one  VG  bar. 
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Figure  25.  (a)  Nusselt  number  versus  mesh  number;  (^)f^Re  number  versus  mesh  number 
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Figure  26.  Mesh  of  the  one  bar  array  domain. 


3.1.6  Fluid  Properties 

The  water  properties  used  in  the  simulations  are  tabulated  in  table  1 .  The  eonstant  properties  were  used  in 
the  eomputational  simulations. 

Table  1  Properties  of  water  used  in  simulations 


Density 

Viseosity 

Heat  eonduetivity 

Speeifie  heat 

998.2  (kg/m^) 

0.001003  (kg/m-s) 

0.6  (w/m-k) 

4182  (j/kg-k) 

3.1.7  Model  Validation 

The  present  numerieal  model  has  been  validated  by  eomparing  the  Nusselt  number  and  pressure  drop 
with  the  experimental  results  available  in  the  literature  [29,  30].  The  eomparison  of  Nusselt  number 
between  numerieal  solutions  and  experimental  results  provided  by  Lee  et  al  [29]  is  shown  in  Figure  27 
using  the  three-dimensional  numerieal  model  at  steady  state.  Figure  28  presents  the  eomparison  of  the 
numerieal  solutions  with  experimental  data  of  pressure  drop  in  mieroehannels  by  Qu  and  Mudawar  [30]. 
The  numerieal  results  for  Nusselt  number  and  pressure  drop  are  in  good  agreement  with  the  experimental 
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results  at  low  Reynolds  number.  However,  when  Reynolds  number  is  larger  than  1500,  the  numerieal 
results  of  Nusselt  number  and  pressure  drop  are  obviously  lower  than  experimental  results.  The  reason 
for  this  is  due  to  the  faet  that  the  laminar  model  was  used  to  simulate  the  fluid  flow  and  heat  transfer 
whereas  the  flow  may  beeome  turbulent  at  high  Reynolds  number  for  mieroehannels. 

To  get  more  eonfidenee  in  the  present  simulation,  the  numerieal  model  has  also  been  tested  for  the  flow 
past  a  eireular  eylinder  plaeed  in  an  otherwise  undisturbed  uniform  eross-flow.  The  Strouhal  number  in 
the  near  wake  of  the  eireular  eylinder  at  Reynolds  number  of  100  was  found  to  be  0.173  well  matehing 
with  the  Srouhal  number  of  0.165  from  published  experimental  results  [32]. 


Figure  27.  Comparison  of  Nusselt  number  between  numerical  results  and  experimental  data 
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Figure  28.  Comparison  of  pressure  drop  between  numerical  results  and  experimental  data 


3.2  Numerical  results  and  discussion 

The  simulations  have  been  earried  out  at  several  different  Reynolds  numbers  eorresponding  to  different 
VG-bar  arrays. 


Numerieal  data  are  ealeulated  using  hydraulie  diameter  defined  as: 


The  Reynolds  number,  Re 


M 


The  Nusselt  number,  Nu 


(18) 


(19) 


(20) 
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And  the  friction  factor,/ 


f  APrf, 

■'  1  2, 

— L 


Where  is  the  axial  pressure  drop. 


(21) 


The  Strouhal  number,  St 


St  = 


V 


Where,  N  is  the  frequeney  of  veloeity  oseillations  and  V  is  the  veloeity. 


(22) 


The  averaged  frietion  faetor  and  Nusselt  number  have  been  tabulated  and  shown  in  Table  2.  The 
tabulated  frietion  faetors  and  Nusselt  numbers  have  been  normalized  with  respeet  to  the  eorresponding 
values  in  the  flow  of  the  ehannel  without  bars  (later  ealled  smooth  ehannel).  From  the  table  one  ean  see 
that  the  frietion  faetors  for  the  mieroehannel  with  bars  are  over  3  times  larger  than  those  of  the 
mieroehannels  without  bars  and  the  Nusselt  numbers  with  bars  are  signifieantly  enhaneed  by  3~5  times  at 
Reynolds  number  of  300  and  400. 

The  effeet  of  the  number  of  VG  bar(s)  on  the  normalized  frietion  faetors  and  Nusselt  number  ean  be  more 
elearly  shown  in  Figure  29  (a)  and  (b).  The  normalized  Nusselt  number  inereased  with  inereasing  number 
of  VG  bars  at  Reynolds  number  of  300.  However,  the  normalized  Nusselt  number  deereased  for  the 
ehannel  with  the  number  of  bars  varied  from  2  to  4  at  Reynolds  number  of  400.  The  normalized  frietion 
faetor  inereased  with  the  number  of  bars  for  both  Reynolds  number  of  300  and  400.  The  normalized 
frietion  faetor  inereased  more  when  the  number  of  bars  ehanged  from  1  to  2  than  when  the  number  of 
bars  ehanged  from  2  to  4. 
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Table  2.  The  averaged  friction  factors  and  Nusselt  number  in  single  microchannel  with  different 
VG-bar  arrays  and  various  Reynolds  number. 


Bar(s) 

Dh  (urn) 

Re 

fbars/f 

NUbars/Nu 

1 

800 

300 

3.03 

3.04 

1 

800 

400 

3.59 

3.80 

2 

800 

200 

3.62 

2.93 

2 

800 

300 

4.71 

3.74 

2 

800 

400 

5.78 

4.88 

4 

800 

300 

5.23 

4.38 

4 

800 

400 

6.01 

4.74 

(a) 


(b) 


Figure  29.  (a)  Non-dimensional  Nusselt  number  for  the  microchannel  with  different  VG-bar  arrays 


at  Reynolds  number  of  300  and  400;  (b)  Non-dimensional  friction  factor  for  the  microchannel  with 
different  VG-bar  arrays  at  Reynolds  number  of  300  and  400; 


35 


Contract  Number:  FA8650-04-C-2405 


Table  3  presents  the  frequeney  of  veloeity  oseillation  and  Strouhal  number  for  various  Reynolds  numbers 
and  bar-arrays.  The  frequeney  inereased  with  the  inerease  of  Reynolds  number.  The  frequeney  with  4 
bars  was  larger  than  the  frequeney  with  2  bars,  whieh  was  larger  than  that  of  1-bar  array.  However,  the 
Strouhal  numbers  for  different  eases  were  very  elose  only  ranging  from  0.99  to  1.1. 

Table  3  The  Frequency  and  Strouhal  number  at  different  Reynolds  number  and  various  VG-bar 
arrays. 


Bar(s) 

Dh  (|im) 

Re 

N  (frequeney) 

St  (Strouhal) 

1 

800 

300 

0 

0 

1 

800 

400 

653.59 

1.05 

1 

800 

600 

925.93 

0.99 

2 

800 

200 

344.83 

1.1 

2 

800 

300 

476.19 

1.02 

2 

800 

400 

ehaotie 

N/A 

2 

800 

500 

ehaotie 

N/A 

4 

800 

300 

483.09 

1.03 

4 

800 

400 

ehaotie 

N/A 

Figures  30  (a),  (b),  and  (e)  show  the  time  traees  of  veloeity  at  Re  =  300,  400,  and  500  respeetively  for 
upstream  and  downstream  points  shown  in  Figure  23  for  the  ehannel  with  2  bars.  The  veloeity  was 
periodieally  ehanging  with  time.  The  unsteady  and  oseillating  behavior  of  flow  was  found.  The  frequeney 
of  veloeity  ean  be  easily  ealeulated  at  Re  =  300  for  the  ehannel  with  2-bar  array.  But  the  flow  was 
ehaotie  and  the  frequeney  was  uneertain  at  Reynolds  number  of  500.  From  these  figures,  we  eonelude 
that  the  veloeity  at  the  upstream  point  is  always  higher  than  that  at  downstream  point,  but  they  have  the 
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same  frequencies  of  oscillation.  The  oscillating  amplitude  of  velocity  of  higher  Reynolds  number  is  larger 
than  that  of  lower  Reynolds  number. 

Figure  3 1  and  Figure  32  show  the  time  traces  of  wall  temperature  and  Nusselt  number  respectively  for  the 
central  point  of  bottom  wall  at  Re  =  300  and  400  with  2-bar  array.  The  oscillating  behavior  of  wall 
temperature  as  well  as  Nusselt  number  with  time  was  found.  The  wall  temperature  increased  and  the 
Nusselt  number  decreased  with  time,  however,  after  four  seconds,  the  averaged  temperature  and  Nusselt 
number  became  constant.  The  averaged  Nusselt  numbers  were  about  10.1  and  13.2  for  Reynolds  number 
of  300  and  400  respectively.  Considering  that  the  Nusselt  number  was  only  about  2.7  for  the  channel 
without  VG  bar,  the  heat  transfer  in  microchannel  with  VG  bars  was  significantly  enhanced  due  to  the 
unsteady  and  oscillating  flows  generated  by  the  VG  bars. 


(a) 
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Figure  33  and  Figure  34  respeetively  show  the  eontours  of  streamwise  and  transverse  eomponent  of 
veloeity  with  different  bar-arrays  at  Reynolds  number  of  400.  The  strongly  unsteady  and  oseillating  flows 
were  found  for  the  ehannels  with  various  bar-arrays.  Results  indieated  that  the  boundary  layers  were 
interrupted  and  beeame  thinner  near  the  VG  bars,  at  the  same  time  the  large-seale  unsteadiness  was  found 
inside  the  mieroehannels.  The  strengths  of  vortieity  with  1-bar  and  2-bar  arrays  were  weaker  than  that  of 
the  ehannel  with  4  bars. 

The  eontours  of  transverse,  spanwise,  and  streamwise  eomponents  of  vortieity  are  shown  in  Figure  35  to 
Figure  37.  It  ean  be  seen  that  the  flow  vortieity  upstream  of  the  bar  is  eonsiderably  stronger  than  the 
vortieity  downstream  of  the  bar.  The  unsteady  natures  of  the  flow  were  indieated  from  the  pattern  of 
vortieity  distribution.  Figure  35  shows  the  eontours  of  spanwise  eomponent  of  vortieity.  It  is  obvious  to 
find  the  formation  of  separating  shear  layers  at  the  leading  edges,  but  the  shear  layers  do  not  roll  up  to 
form  large-seale  vertieal  strueture  behind  the  bars.  The  transverse  eomponent  of  vortieity  is  presented  in 
Figure  36,  it  also  shows  that  the  upstream  vortieity  is  mueh  stronger  than  downstream  vortieiy.  Although 
the  vortieity  in  transverse  direetion  was  weak,  it  may  have  played  an  important  role  in  enhaneing  heat 
transfer  due  to  eonsiderable  inerease  in  Nusselts  number  shown  in  Table  2.  From  Figures  35  and  36,  the 
ehannel  with  4-bar  array  has  the  strongest  vortieity.  The  behavior  of  streamwise  vortieity  ean  be  found 
form  Figure  37,  it  also  shows  the  ehannel  with  more  bars  has  stronger  vortieity. 
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Figure  33.  Contours  of  spanwise  component  of  velocity  with  different  bar-arrays  {Re  =  400, 
Dh  =  0.8mm,  time  =  0.8s) 
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Figure  34.  Contours  of  transverse  component  of  velocity  with  different  bar-arrays  {Re  = 
400,  Dh  =  0.8mm,  time  =  0.8s) 
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Figure  35.  Contours  of  spanwise  component  of  vorticity  with  different  bar-arrays  {Re  = 
400,  Dh  =  0.8mm,  time  =  0.8s) 
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Figure  36.  Contours  of  transverse  component  of  vorticity  with  different  bar-arrays  (Re  = 
300,  Dh  =  0.8mm,  time  =  0.18s) 
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Figure  37.  Contours  of  streamwise  component  of  vorticity  at  mid  plane  with  different  bar- 
arrays  (Re  =  400,  time  =  0.  8s):  a)  1-bar  array;  b)  2-bar  array;  c)  3-bar  array. 


Figures  38  and  39  present  the  temperature  eontours  of  bottom  wall  with  1-bar  and  4-bar  arrays 
respeetively.  Figure  41  and  Figure  42  show  the  three-dimension  temperature  eontours  at  the  bottom  wall 
and  the  surfaees  of  bars.  Although  the  downstream  and  upstream  veloeities  of  the  ehannel  with  4  bars 
were  weak  eompared  with  other  ehannels,  the  temperature  of  bottom  wall  of  the  ehannel  with  4  bars  was 
more  uniform  and  lower  than  that  with  one  bar  due  to  high  heat  transfer  resulted  from  its  strong  vortieity 
and  mixing.  This  means  the  ehannel  with  more  bars  has  enhaneed  the  heat  transfer  better.  The  ehannel 
with  4-bar  array  ean  sustain  larger  heating  flux  at  bottom  wall.  Figure  40  shows  the  temperature  eontour 
of  bottom  wall  with  2-bar  array  but  the  heat  flux  was  set  to  200,000  w/m^,  while  the  heat  flux  was  set  to 
50000  w/m^  for  the  ehannel  with  1  and  4  bars.  With  this  heat  flux  the  temperature  of  bottom  wall  was 
signifleantly  inereased  and  more  uneven,  espeeially  at  the  edges  of  the  bottom  wall. 


It  is  interesting  to  find  that  the  temperature  at  the  region  near  the  bars  was  lower  than  that  at  the  other 
area  of  the  heating  wall.  The  VG  bars  improved  the  heat  transfer  in  mieroehannels  due  to  the  unsteady 
and  oseillating  fluid  flow. 
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Figure  38.  The  temperature  contour  at  bottom  plane  {Re  =  400, 1  bar,  time  =  0.8  sec,  q  = 
50000  w/m^) 
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Figure  39.  The  temperature  contour  at  bottom  plane  (Re  =  400, 4  bars,  time  =  0.8  sec,  q  = 
50000  w/m^) 
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Figure  40.  The  temperature  contour  at  bottom  plane  {Re  =  400,  2  bars,  time  =  0.8  sec,  q  = 
200000  w/m^) 
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Figure  41.  The  3-D  temperature  contour  at  bottom  plane  and  surfaces  of  bar  {Re  =  400, 1 
bar,  time  =  0.8  sec,  q  =  50000  w/m^) 


Figure  42.  The  3-D  temperature  contour  at  bottom  plane  and  surfaces  of  bar  {Re  =  400,  4 
bars,  time  =  0.8  sec,  q  =  50000  w/m^) 
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3.3  Conclusions  from  the  numerical  study 

MicroChannel  with  VG  bars  can  provide  high  heat  transfer  surface  area  per  unit  volume  and  alter  the  fluid 
dynamics  to  enhance  mixing.  Numerical  study  of  fluid  flow  and  heat  transfer  in  the  microchannel  with 
VG  bars  resulted  in  following  findings: 

•  The  flow  in  microchannel  was  found  to  be  unsteady  and  oscillating  at  Reynolds  number  equal  to 
or  larger  than  300  with  2  bars  and  4  bars.  For  the  microchannel  with  1  bar  the  flow  was  found  to 
be  naturally  steady  at  Reynolds  number  of  300.  At  Reynolds  number  of  500  with  2  bars,  the  flow 
in  the  microchannels  could  become  chaotic. 

•  For  most  of  the  microchannels  with  bars,  the  friction  factor  was  more  than  3  times  larger  than 
that  for  the  microchannels  without  bars  at  the  same  Reynolds  number. 

•  The  larger  the  Reynolds  number  and  the  more  the  number  of  VG  bars,  the  stronger  the 
unsteadiness  and  oscillations  of  the  flows 

•  With  increasing  Reynolds  number  the  averaged  temperature  of  heated  wall  was  decreased 

•  The  oscillating  behavior  of  velocity,  averaged  temperature  on  the  bottom  wall,  and  Nusselt 
number  was  found  at  Reynolds  number  of  300  or  larger  for  the  channel  with  2  or  4  bars,  and  at 
Reynolds  number  of  400  or  larger  for  the  channel  with  1  bar. 
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